###########   Cervejas

### 1-) Representao de cervejas em grfico de barras e dotchart.
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
table(cevejas)
barplot(table(cervejas))

### 2-) O dotchart possui uma maior razo dado/tinta.

###########  Rivers

### 1-)
prop.rios.peq = (length(rivers[rivers<mean(rivers)]))/length(rivers)
### 2-) 3 quantil = 680
summary(rivers)
  Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  135.0   310.0   425.0   591.2   680.0  3710.0 
### 3-)
mean(rivers, trim=0.25)
449.9155
### Mdia = 591.2, mdia truncada(0.25)= 449.9155  diferem, a ultima  mais prxima da mediana (425.0), pois valores extremos esto puxando a mdia para cima.

########### CAIXETA
### 1-) 
hist(caixeta$fuste)
### 2-)
retiro = caixeta[caixeta$local=="retiro",]
jureia = caixeta[caixeta$local=="jureia",]
chauas = caixeta[caixeta$local=="chauas",]
hist(retiro$h, )
hist (jureia$h, histogram(~ h | local, data=caixeta))
hist (chauas$h)
### 3-) 
summary(jureia$h)
summary(retiro$h)
summary(chauas$h)
Pelos histogramas, observa-se que as classes de alturas mais frequentes esto entre 50-150, nas trs localidades. Chauas possui rvores entre 450 - 500, 
que no ocorrem nas outras reas. Mesmo assim, a mdia de altura das rvores de jureia  maior e a de retiro a menor das trs.

########### E.grandis
eucalipto <- read.table("egrandis.csv", header=TRUE, sep= ";")
str (eucalipto)
summary(eucalipto)
### 1-)
boxplot(dap~regiao, eucalipto)
boxplot(dap~rotacao, eucalipto)
### 2-)
quantile(eucalipto$h, probs= seq(from=0,to=1,by=0.1))
qqnorm(eucalipto$h)
qqline(eucalipto$h)

############ Caixeta 
###Isolando Caixetas
s.caixetas <- caixeta[caixeta$especie=="Tabebuia cassinoides",]
###Transformando CAP em DAP
s.caixeta$DAP = s.caixeta$cap/pi
xyplot(dap~h|local,data=tabebuia)
###
par(mfrow=c(2,2))
plot(DAP~h, data=s.caixeta,subset=local=="retiro", main="T. cassinoides de Retiro",xlab="Altura",ylab="DAP")
plot(DAP~h, data=s.caixeta,subset=local=="chauas", main="T. cassinoides de Chauas",xlab="Altura",ylab="DAP")
plot(DAP~h, data=s.caixeta,subset=local=="jureia", main="T. cassinoides de Juria",xlab="Altura",ylab="DAP")
### Linearidade funo Smooth
par(mfrow=c(2,2))
scatter.smooth(s.caixeta$DAP[s.caixeta$local=="retiro"]~s.caixeta$h[s.caixeta$local=="retiro"], main="T. cassinoides de Retiro",xlab="Altura",ylab="DAP")
scatter.smooth(s.caixeta$DAP[s.caixeta$local=="chauas"]~s.caixeta$h[s.caixeta$local=="chauas"], main="T. cassinoides de Chauas",xlab="Altura",ylab="DAP")
scatter.smooth(s.caixeta$DAP[s.caixeta$local=="jureia"]~s.caixeta$h[s.caixeta$local=="jureia"], main="T. cassinoides de Jureia",xlab="Altura",ylab="DAP")
### Pacote Lattice
library(lattice)
xyplot(dap~h|local,data=tabebuia)
